Magnetic edge states in spin triplet superconductors 
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We show that a spontaneous magnetic moment may appear at the edge of a spin-triplet supercon- 
ductor if the system allows for pairing in a subdominant channel. To unveil the microscopic mecha- 
nism behind such effect we combine numerical solution of the Bogoliubov-De Gennes equations for 
a tight-binding model with nearest-neighbor attraction, and the symmetry based Ginzburg-Landau 
approach. We find that a potential barrier modulating the electronic density near the edge of the 
system leads to a non-unitary superconducting state close to the boundary where spin-singlet pairing 
coexists with the dominant triplet superconducting order. We demonstrate that the spin polariza- 
tion at the edge appears due to the inhomogeneity of the non-unitary state and is manifested via 
lifting of the spin-degeneracy of the Andreev bound-states. 

PACS numbers: 74.20.Rp,74.25.Dw,74.70.Pq,71.10.Li 



Introduction. Recognition that the surface states in 
correlated materials reflect the nature of the interactions 
and orders existing in the bulk has led to a very sig- 
nificant research effort aimed at the understanding, and 
potential control, of these electronic states Sys- 
tems where the bulk of the material is gapped but the 
boundary supports gapless modes are especially interest- 
ing since under these conditions the surface states are ro- 
bust, and may be topologically protected, i.e. their exis- 
tence relies on the global symmetries of the bulk state and 
does not depend on the details of the surface scattering 
and other sample-dependent parameters Q. The bulk 
gap may be due to the band structure, or, in a metal, may 
arise at low temperatures from electron-electron interac- 
tion, as in superconductors {3]. Simple band insulators 
or conventional superconductors do not support robust 
low-energy states at the boundary, but it is the study 
of their counterparts with the bulk that is topologically 
non-trivial, and hence the bulk-boundary correspondence 
theorem dictates the existence of the surface states, that 
has been a focus of much recent attention 0, Q, H| • 

One of the best candidates for the topological super- 
conductivity is Sr 2 Ru04, where the emergent consensus 
indicates triplet chiral pairing, with time-reversal sym- 
metry broken by the orbital degrees of freedom [f|. In 
this material the topologically protected edge states have 
been predicted, and their signatures were recently found 
in tunnelling spectroscopy [6| . The quasiparticles reflect- 
ing off the sample boundary experience the sign change of 
the superconducting order parameter along their trajec- 
tory, which gives rise to so-called Andreev bound states 
(ABS) near the surface. Emergence of ABS has been well 
investigated in high-T c cuprates and a number of other 
unconventional superconductors 0, 0, H| ■ 

In this Letter we investigate the nature of the Andreev 
bound states at the surface of spin triplet superconduc- 
tors. We perform a microscopic self-consistent calcu- 



lation, including the realistic surface barrier of a finite 
width and height and the possibility of pairing in one or 
more subdominant channels 0- We find that a) a sub- 
dominant in-phase s-wave superconducting order exists 
near the edge of the sample; b) the in-phase s-wave com- 
ponent gives a non-unitary superconducting state at the 
boundary; c) MS M result, the ABS is spin-polarized, lead- 
ing to a finite surface magnetization. We analyze the con- 
ditions for the existence of the magnetic surface states, 
and we investigate their spectrum. These numerical re- 
sults are confirmed by carrying out a Ginzburg-Landau 
expansion of the free energy. Our work strongly suggests 
that triplet superconductors can be used in spin-active 
heterostructures. 

Model and formalism. We consider a two-dimensional 
superconductor in a slab geometry with two parallel in- 
terfaces separating it from the vacuum. Denoting by x 
and y the directions perpendicular and parallel to the in- 
terfaces, respectively, we assume that the system is uni- 
form along the y axis, so that the translational symme- 
try is broken only in the x direction. The Hamiltonian 
is then defined on a square lattice of size L x L (we set 
the lattice constant to unity), with periodic boundary 
conditions along y, 

h = -t ( c L c i<y + h - c -) - l- l J2 n ^ 

~J2 V ( n n n il + n U n rt) +^2 U {ix) n ™ ■ (!) 
<ij> i 

Here the lattice sites are labelled by i = {i x ,i y ), with i x 
and i y integers between and L, denote nearest- 

neighbor sites, and [i is the chemical potential. The 
nearest-neighbor attractive interaction — V (V > 0) is 
effective in both singlet and triplet pairing channels. All 
the energies are in units of the hopping parameter t. The 
two parallel edges of the system are located at i x = 
and i x = L, and we introduce a site-dependent potential 
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U(i x ) to model the interface barrier. To investigate the 
model of Eq. ([I]) we decouple the interaction term in the 
Hartree-Fock approximation by introducing the pairing 
amplitude on a bond, Ay = (ci^Cj^), so that Vn-^n^ ~ 
(Ayc]j_4 t + A y c it Cj4, - |A ;j | 2 ). These pairing ampli- 
tudes on each bond yield the spin singlet (S) and triplet 
(T) components, A S < T = (Ay ± AjO/2, that define the 
superconducting order parameters with s- or p-wave sym- 
metry, i.e. A s (i) = (Af i+ - + A?j_ s + A? i+ - + Af,_ y )/4 

and A /M B )(i) = ( A U+x(y) - ^li-x(y))/ 2 ' Which arC th(3n 

determined self-consistently [9|. Singlet d-wave super- 
conductivity is possible but does not appear in the pa- 
rameter range where we work Q. In the bulk (U = 0) 
the most favorable pairing state for this model depends 
on the electron density, n, and, in particular, the chiral 
Px + ipy order is stabilized in the region between half- 
filling, (i ~ 0, and high (low) density (|/i| ~ 2.5) [Toj j . 
Hence, we choose the chemical potential to be ~ 1.8, 
in the window of stability, so that in the absence of U(i x ) 
the filling is n w 0.4. All the numerical results below have 
been obtained for a pairing interaction V = 2.5 and a sys- 
tem size L = 80, and we checked that greater values of L 
leave the results qualitatively unchanged. Many studies 
of the surface states assume a sharp step-like potential at 
the surface [Uj]: we find qualitative differences between 
the results obtained using such assumptions and the be- 
havior of the system using a realistic surface barrier. One 
very important distinction is that a finite-width barrier 
changes the electron density near the boundary, thereby 
enabling the emergence of superconducting components 
competing with the dominant triplet one. 

Numerical results. FigfT] shows representative results 
for the electron density and spin polarization, as well 
as the evolution of the superconducting order parame- 
ters for different strength of the surface potential. Here 
we assumed a rectangular potential barrier of height U 
(in units of the hopping t) near the left edge of the sys- 
tem, < i x < 8. Finite values of U lead to the deple- 
tion of the electron density near the edge, Figflja). For 
U = 0, FigQTc), we find the expected result: the inter- 
face is pairbreaking for the p x component of the order 
parameter, while the p y component remains essentially 
constant. However, as U exceeds a critical value, here 
found to be U c ~ 0.19, a subdominant s-wave compo- 
nent of the order parameter appears, Figjljd)-(f), and 
there is a substantial region of coexistence of the super- 
conducting order parameters of different parity. Note 
that mixed parity order parameters are generally allowed 
here since the presence of the barrier breaks the inversion 
symmetry of the lattice. A remarkable feature is that the 
emergence of the mixed-parity phase is accompanied by 
the appearance of a finite spin polarization in that same 
region, FigQJb). As the barrier becomes higher, fewer 
carriers remain in the boundary layer, and the magneti- 
zation gradually decreases. 
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FIG. 1. (color online) Evolution of the electron density and 
the order parameters with the surface barrier. The s-wave 
refers to nearest-neighbor pairing, a-b) spin-up electron den- 
sity and spin polarization for the several values of the barrier 
height U, c)-f) superconducting order parameters for different 
heights of the potential barrier U extending from i x = to 
ix = 8. 

Analysis of the energy spectrum E n (k y ) (due to trans- 
lational invariance along the interface k y is a good quan- 
tum number) obtained from the numerical solution of 
the Bogoliubov-De Gennes equations confirms that the 
local magnetization is due to the gapless edge modes 
propagating in one direction along the system bound- 
ary. This is clear from Fig. [2J where the two originally 
spin-degenerate chiral edge states associated with the 
left boundary become split once the barrier potential ex- 
ceeds U c ~ 0.19. This splitting occurs concomitantly 
with the appearance of a spin polarization at that edge 
(of course, the countcrpropagating edge mode associated 
with the right edge is not affected) . As the barrier height 
is increased the splitting between the spin-up and spin- 
down modes at the edge decreases, and eventually leads 
to lower and lower values of the magnetization, as the 
boundary region becomes depleted. We also notice that 
extra midgap states appear close the bottom and the top 
of the gap edge with a structure that evolves from an 
asymmetric profile in the case of an extended region of 
singlet-triplet coexistence at the surface, i.e. for U = 0.2 
and 1.0, to a more symmetric one around k y = for 
U = 2.2. 
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The implication of this result is that a spin accumu- 
lation may occur in a triplet superconductor without the 
proximity coupling to an exotic system. The characteris- 
tic length scale over which the magnetization appears 
in our case is comparable to the width of the region 
where the density gradient is present: this is because 
for the present values of parameters the superconduct- 
ing coherence length is comparable to that width. In 
general, if the screening of the surface potential occurs 
on a much shorter scale than the coherence length, the 
surface barrier nucleates the subdominant s-wave com- 
ponent of the pairing amplitude, and the length scale of 
the coexistence is set by the coherence length. Note that 
this situation is quite different from the case discussed 
m RefsEMl where spin polarization appears at the in- 
terface separating a triplet and a singlet superconductor 
only for a non-trivial phase difference. In fact, as we show 
below, the requirements for the spin accumulation are 
very different for the two geometries. The physical dis- 
tinction is that the studies of superconductor-insulator- 
superconductor (S-I-S) junctions within standard BTK- 
like techniques use the scattering formulation of the 
boundary problem between two semi-infinite materials, 
whereas here the polarization is obtained from a fully self- 
consistent solution of the Bogoliubov-de Gennes equa- 
tions in the presence of a realistic boundary potential. 

Symmetry analysis. To elucidate the origin of the spin 
polarization we analyse the problem from the symmetry 
perspective. The main insight from the numerical results 
is that the singlet and the triplet components of the order 
parameter coexist over a finite length scale near the edge 
due to the extended nature of the boundary potential, 
and because of the absence of the inversion symmetry 
near the interface. We first note that in this situation the 
combined order parameter is non-unitary. In the 2x2 spin 
space the order parameter is commonly written in the 
form [lH Afc = i [(d k ■ <r) + ip k ] °~y, where o~i arc Pauli 
matrices, and dk and ip k are the three triplet and the 
singlet pairing amplitudes respectively. In that case the 
gauge-invariant product is given by AA^ = \d k \ 2 o~o + 
iV'fcl 2 'Ik ■ rr with q k = ijj k d* k + i\) k d k + i[d k x d" k . 
In our case the bulk triplet superconductor is unitary, 
d\\z and therefore d k x d k = 0. However, in the coexis- 
tence region, q k = 2Rc [ipkd k ] ^ since, as Fig. [T]shows, 
the self-consistent solution yields the in-phase singlet and 
triplet p x components. 

Nonvanishing q k vector suggests, but by itself docs not 
necessarily require, a finite spin polarization, since its av- 
erage over the direction of the momenta fc may still van- 
ish (this averaging is difficult to carry out analytically 
in a system without translation invariance such as our 
boundary problem). Therefore we consider a Ginzburg- 
Landau (GL) expansion of the magnetic contribution to 
the free energy density f m in the region of coexistence. 
We choose the direction of the magnetization m paral- 
lel to d(r) = z (rj x (r) + ii] y (r)), along the £-axis, as rc- 




o 

<l 



M 

m 



- U=0.2 / r """^\ s 






Sink, ,- 



1.5 
1 

0.5 


-0.5 
-1 
-1.5 
1.5 
1 

0.5 


-0.5 
-1 
-1.5 
1.5 
1 

0.5 


-0.5 
-1 
-1.5 
1.5 
1 

0.5 



-0.5 
-1 
-1.5 



FIG. 2. (color online) Energy spectrum of the solutions of the 
Bogoliubov-de Gennes equations scaled to the superconduct- 
ing gap Ao of the homogeneous spin triplet superconductor. 
Upper panel: pure chiral phase. Two dispersing branches 
crossing E = correspond to the modes along the left and 
the right edges of the sample. Lower panels: mixed-parity 
phases. The right edge mode remains unchanged while the 
modes along the left edge become spin-split. Extra midgap 
states which are spin-polarized appear at the bottom (occu- 
pied) and the top edge (unoccupied) of the gap reflecting the 
occurrence of a mixed-parity coexisting region. The latter 
and the spin-split chiral states both contribute to the overall 
spin accumulation. 



quired by the spin rotation invariance, and assume that 
the pairing amplitudes depend only on the coordinate 
x normal to the boundary. Then the GL expansion al- 
lows for the terms linear in m and the gradient of the 
Pa-component of the triplet pairing, 




fm = m 2 + {d x m) 2 



f am^/3d x [r) x ip* + n x 4> 



(2) 



Since the presence of the linear in m term means that 
for any values of the coefficients the minimum of the free 
energy is at finite to, it immediately follows that in the 
region of the coexistence a finite magnetization appears 
unless the singlet and the triplet p x components are out 
of phase and the product j] x ij} is purely imaginary. This 
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observation emphatically brings forth the distinction be- 
tween our results and those for the S-I-S junction [l2- 14[. 
where the spin accumulation only occurs if the two pair- 
ing components are out of phase, the exact opposite of 
the result we find. 

Only the p^-component of the triplet appears in the 
GL expansion above; in principle the term d y r\ y is also 
allowed by symmetry , but does not appear under the 
assumption of the translational invariance along the in- 
terface. It follows that the time-reversal symmetry break- 
ing by the bulk chiral triplet state is not at the origin of 
the magnetization of the Andreev bound states: the same 
result would be achieved for purely real p x bulk triplet su- 
perconductivity, while for the imaginary p x bulk pairing 
with real subdominant s-wavc pairing near the interface 
no magnetization appears. To explore this connection, 
we consider a situation when in proximity of the left edge 
(0 < i x < i Xl with i x = 10) a purely local s-wave poten- 
tial is the only source of pairing, with V still effective 
in the remaining part of the system, and use as input a 
real s-wave order parameter for < i x < i x and a purely 
imaginary p^-wave one for i x > i x . Working at the same 
electron density as before (with U = 0), and without the 
self-consistent iterative procedure (which would unavoid- 
ably lead to a mixing of real and imaginary components), 
we see from the panels a) and c) of Fig[3]that, in spite 
of the parity mixing occurring around i x = i x , no appre- 
ciable spin polarization is observed. On the other hand, 
when in the same configuration one uses for i x > i x a 
purely real p x input order parameter to generate a non- 
unitary mixed phase around i x = i Xl a significant mag- 
netization clearly develops in the mixing region. This is 
exhibited in the panels b) and d) of Fig. [3J and supports 
our conclusion based on symmetry arguments. 

Discussion. We show that Andreev bound states near 
a boundary of a triplet superconductor can be spin- 
polarized. The origin of the spin polarization is in the 
emergence of the coexistence regime of the triplet and the 
subdominant singlet pairing components near the inter- 
face. Our numerical results demonstrate that the two are 
phase-locked, and both the numerical fully self-consistent 
solution of the Bogoliubov-de Gennes equations and the 
Ginzburg-Landau analysis indicates that magnetization 
inevitably appears when the two order parameters lead 
to a non-unitary configuration and arc spatially vary- 
ing. We find that the symmetry-breaking at the surface 
is unrelated to the chiral nature of the bulk supercon- 
ducting state, and therefore may be expected in a much 
wider class of triplet superconductors. It would also be 
very interesting to check whether similar effects occur at 
the interfaces involving non-centrosymmetric supercon- 
ductors, where the singlet and the triplet components are 
intrinsically mixed in the bulk yielding measurable spin 
effects at the interface 17 -2(|, as well as in the proximity 
structures with topological materials. We leave this for 
future investigations. 
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FIG. 3. (color online) Spatial evolution of singlet and triplet 
order parameters and of the magnetization in the case of uni- 
tary, a) and c), and non-unitary states, b) and d), respec- 
tively, as realized in a thin vertical slab of the system around 
the point (i x = 10). A real spin singlet on-site pairing ampli- 
tude and a purely real or a purely imaginary triplet pairing 
amplitude are assumed in input (i x < 10 and i x > 10, respec- 
tively) . 
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